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We consider a curved chain of nonlinear oscillators and show that the interplay of curvature and 
nonlinearity leads to a symmetry breaking when an asymmetric stationary state becomes energeti- 
cally more favorable than a symmetric stationary state. We show that the energy of localized states 
decreases with increasing curvature, i.e. bending is a trap for nonlinear excitations. A violation of 
the Vakhitov-Kolokolov stability criterion is found in the case where the instability is due to the 
softening of the Peierls internal mode. 



Nanoscale physical systems such as biological macro- 
molecules, nanotubes, microtubules, vesicles, electronic 
and photonic wave-guide structures Jl]-|| are usually ob- 
jects with nontrivial geometry. The gene transcription 
is usually accompanied by a local DNA bending 
the electronic properties of carbon nanotubes drastically 
depend on their chirality JTo[ ] ; T-shape junctions were re- 
cently proposed |p"l|| as nanoscale metal-semiconductor- 
metal contact devices. The geometry manifests itself in 
particular in creation of linear quasi-bound states (see 
e.g. 1 12 13]). On the other hand, the competition be- 



tween nonlinearity and dispersion likewise results into 
localization of energy. Taken together these two localiza- 
tion mechanisms are competitive and one might expect 
that the interplay of geometry and nonlinearity can lead 
to qualitatively new effects. An example of nonlinearity- 
induced change of the geometry of the system is found 
in Ref. [03 where the classical Heisenberg model on a 
two-dimensional manifold was considered. It was shown 
that a periodic topological spin soliton induces a periodic 
pinch of the cylindric manifold. 

In this paper we consider nonlinear Schrodinger (NLS) 
models on a curved one-dimensional manifold and show 
that the competition between curvature and nonlinear- 
ity leads to a qualitatively new effect: curvature-induced 
symmetry breaking for non-topological solitary excita- 
tions. This effect is inherent for both continuum and 
discrete NLS models. We start with the discrete NLS 
equation 



d , 
dt 



(1) 



which describes a system of coupled nonlinear oscilla- 
tors, ip n (t) being the complex amplitude at the site n 
(n = 0, ±1, ±2). In the theory of charge (energy) transfer 
ip n (t) is the wave function of the carriers. The Hamilto- 
nian 



and the number of excitations (quanta) N = ^ n \il> n \ 2 
are conserved quantities. In Eq. ([[]) the nonlinear term 
|'0n| 2 '0n represents a self-interaction of the quasiparticle 
(anharmonicity of vibrations). The excitation transfer 
Jnm in the coupling term depends on the distance in the 
embedding space between the sites n and m: J nm = 
J(\r n - f m \), where the radius-vector r„ = [x n ,y„,z n ) 
characterizes the position of the site n on the mani- 
fold. The sites are assumed to be equidistantly placed on 
plane envelope curve (Fig. |l|). This model is closely re- 
lated to the so-called wormlike incxtensible chain model 
Jl5| , ^6| which is abundantly used in polymer dynamics 
(see e.g. JT7|-[T9[] ) . We consider the properties of nonlin- 
ear excitations in the vicinity of a smooth bend where 
the curve can be modeled by the parabola without loss 
of generality (as done in Fig. [l]): y n = kx^/2 and 
z n = with k being the inverse radius of curvature at 
the bending point. The distance between neighboring 
sites is assumed to be independent on the curvature k: 
\r n +i — r n \ = 1 (inextensibility of the chain). There- 
fore, the excitation dynamics would not depend on the 
curvature of the system in the nearest-neighbor approxi- 
mation, when J nm = J8 n -m.±i- In many physical prob- 
lems, however, the nearest-neighbor approximation is too 
crude. For example, the DNA molecule contains charged 
groups, and therefore the vibration-excitation transfer is 
due to the dipole-dipole interaction decaying with the 
distance r as 1/r 3 . Effective dispersive interaction be- 
tween nonlinear layers in nonlinear dielectric superlat- 
tices is exponential, exp(— ar), with the inverse radius 
a depending on parameters of the lattice (superlattice 
spacing, linear refractive index and so on) pOf . With 
this in mind one can see that in general case the num- 
ber of neighbors available for excitation is larger near the 
bending point than far away from it. This gives birth, 
in particular, to an effective double-well potential U (see 
Fig. [I] and Eq. (|^)) which is basically responsible for the 
effects discussed below. 

We investigated both J nm = Jexp (— a\f n — f m \) and 
Jnm = J\r n — T m | ~ s . But in this paper we discuss in 
detail the first case only. It is known (21j that in the 
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case of straight line (k = 0) such NLS model exhibits 
(for a < 1.7) bistability in the spectrum of nonlinear 
stationary states. To distinguish these effects from the 
finite curvature effects we use a = 2. We study station- 
ary states of the system ip n (t) = <f> n (A) e tAt , where A is 
the nonlinear frequency and <p n {A) is the amplitude in 
the n-th site. 




FIG. 1. Curved chain with sites equidistantly placed on a 
parabola with the curvature k = 0.5 at x — and correspond- 
ing effective double-well potential U n = — 53 Jn,m 
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FIG. 2. The number of excitations N(A) for the straight 
chain (k = ; dot-dashed line) and curved chain (k = 3 ; 
solid, dotted, and dashed lines) for a — 2 and J = 6.4. In 
the inset shapes of the symmetric (s) and asymmetric (as) 
stationary states for N = 3.6 are shown. 

Figure || shows N versus A obtained numerically for 
the straight chain and curved chain. Several new fea- 
tures arise as a consequence of the finite curvature: 

• There is a gap at small A in which no localized 
excitation can exist. This gap originates from the 
existence of a curvature created linear two-hump 
localized mode. The gap increases when the curva- 
ture k increases. 

• There are two branches of stationary states: a 
branch of symmetric localized excitations (s), 
which exists for all values of N, and a branch 
of asymmetric localized excitations (as), which ex- 
ists only for N > Nth(K). The threshold value 
N t h(n) decreases when n decreases and vanishes in 



the limit k — > 0. The symmetric stationary state 
has a two-humped shape (evolved from the linear 
two-hump localized mode) while in the asymmetric 
state the maximum is shifted either to the left or to 
the right from the center of symmetry n — (see 
inset in Fig. ||). 

For the symmetric stationary states N(A) is non- 
monotonic as in Refs. for NLS models on 
the straight chain with the radius of the disper- 
sive interaction above a critical value. Similarly, in 
the case of curved chain there is an interval of N 
in which three different symmetric states coexist 
for each excitation number. Taking into account 
that for the inverse radius a = 2 under investiga- 
tion there is no bistability in the straight chain (see 
dot-dashed line in Fig. H), one may conclude that 
bistability is facilitated in the systems with finite 
curvature. 




FIG. 3. Spectrum of linear excitations around the symmet- 
ric stationary state <p n (A) for a = 2, J = 6.4 and k = 3. There 
are two internal localized modes: Peierls antisymmetric (p) 
and breathing symmetric (b). In the inset the dependence 
N(A) is shown. 

The linear stability of the stationary states is in- 
vestigated (as it is described in Refs. p3| , ^l[ ) letting 
ipn(t) — [4>n(A) + e„exp(ifM)] exp(iAi) and linearizing 
Eq. (Q) with respect to e„. By definition the station- 
ary state </> ra (A) is linearly stable if fl 2 > 0. Figure || 
shows the spectrum of linear excitations around the sym- 
metric localized state. We find that there are two local- 
ized internal modes: Peierls antisymmetric mode (p) and 
breathing symmetric mode (b). There are two intervals 
of instability of the symmetric stationary state: when 
0.99 < A < 1.26 the instability is due to the breathing in- 
ternal mode. In this part of the curve N(A) (dashed line 
in the inset) dN/dA < and the instability agrees with 
the Vakhitov-Kolokolov stability criterion |24|| . When 
0.84 < A < 0.99 the stationary state is unstable with 
respect to the Peierls mode (the corresponding part of 
N(A) is marked by the dotted line). In this latter case 
we observe a violation of the Vakhitov-Kolokolov stability 
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criterion because this instability interval corresponds to 
a part of the curve TV(A) where dN/dA > (see Fig. ||). 




FIG. 4. Energy of symmetric (s) and asymmetric (as) non- 
linear excitations versus the curvature k for a — 2, J = 6.4, 
and TV = 4. 
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FIG. 5. Switching from the symmetric state to the asym- 
metric state at a = 2, J = 6.4, k = 4, and TV = 3.96. 



The energy of the symmetric as well as asymmetric lo- 
calized states is monotonically decreasing function of the 
curvature (see Fig. ||). Thus, one may expect that the 
nonlinear localized excitation may facilitate bending of 
flexible molecular chain. For N > N t h(n), when symmet- 
ric and asymmetric stationary states coexist, the asym- 
metric state is always energetically more favorable. Tak- 
ing into account that in the linear limit the ground state 
of the model is symmetric with respect to the center of 
symmetry n = one can conclude that combined ac- 
tion of the finite curvature and nonlinearity provides the 
symmetry breaking in the system. A typical evolution 
of the curved NLS model in the case when the symmet- 
ric stationary state becomes unstable is shown in Fig. 
§. One can see that the two-humped initial symmetric 
state evolves into an asymmetric state with an excited 



breathing internal mode. 

Apart from the smooth bending we studied the case 
of the edge: y n = \x n \ tan(0) . Qualitatively the same 
results were obtained for this case, too. 
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FIG. 6. The dependence TV (A) for symmetric (s) and asym- 
metric (as) nonlinear excitations in the two-impurity minimal 
model @ with e = 0.7 and a = 1. 

To clarify the physical meaning of obtained results we 
pass to the continuum limit of Eq. (||) where the Hamil- 
tonian H can be represented as 



H 



1 



L 2 W 



9-0 



d£ 



■ U(£)\ip\ 



1 



m 4 \d£. (3) 



Here I is the arclength (the continuum analogue of n), 
L 2 {£) = J^i' - if J u ,d£ is the square of the disper- 
sive length of the excitation, and U(£) = — Je^'di' 
is the energy shift due to the coherent interaction. The 
inverse length L _1 (£) and attractive potential U(£) are 
enerally double-well functions. At the bottom of Fig. 

shown. In the vicinity of the bend there are two sym- 
metric sections of the chain in which the number of 
neighbors entering into the range of the dispersive inter- 
action, 1/q , is higher than anywhere else on the chain. 
In this way a double- well potential structure results. The 
positions and depths of the wells are determined by the 
range of the dispersive interaction and the peculiarities 
of the bending geometry. E.g., in the case of the edge, 
Vn = \ x n\ tan(0) , we get for small 9 

L 2 (£) » £ (4 + ((a 2 \l\ + l) 2 + l)( e ^\i\cos(26) _ e -a\e\^ 



-a\i\ cos(20) 



-a\t 



and U{£) « (2 
It is worth noting that the Hamiltonian similar to Eq. 
(|) with U{1) = was introduced in Ref. @ to study 
the effects of the DNA bending on breather trapping. 
The step function dependence of the dispersive length 
L(£) was there postulated. In all cases bending manifests 
itself as a trap for nonlinear excitations. To obtain an 
analytical solution of the problem one can approximate 
L(£) by a constant and U(£) by two ^-functions. Then 
we arrive at the following minimal model 
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i ~^t + ~di + e{6{e ~ a) + 5{e + a)) ^ + ^ = ' (4) 

where the parameter e is the depth of the well. Intro- 
ducing into Eq. (|4|) the stationary state ansatz ip(i,t) — 
4>(l) exp(iAi) we find that the corresponding nonlinear 
eigenvalue problem can be easily reduced to an algebraic 
one. The results of the analytical consideration are shown 
in Fig. ^. The stationary state (s) is unique and symmet- 
ric for small number of quanta N . But when the num- 
ber of quanta exceeds some critical value there are two 
stationary states: symmetric (s) and asymmetric (as), 
with the asymmetric state being energetically favorable. 
Thus, similarly to the original model, the minimal model 
demonstrates a symmetry breaking effect (localization in 
one of the wells). 

We conclude that in the NLS model the interplay of 
curvature and nonlinearity leads to the symmetry break- 
ing when the asymmetric stationary state becomes en- 
ergetically favorable than the symmetric one. We have 
found that the energy of localized states decreases with 
the increasing of curvature. One may expect that the 
nonlinear localized excitation may facilitate bending of 
a flexible molecular chain. In all cases bending man- 
ifests itself as a trap for nonlinear excitations. For the 
curved NLS model in the case of the physically important 
dipole-dipole dispersive interaction J nm = \r n — r m \~ 3 
qualitatively the same properties are found. Using as 
an example the Davydov-Scott model ^6|,^7j with the 
dipole-dipole coupling between sites (2{| for the param- 
eters which characterize the motion of the amide-I exci- 
tation in proteins (they correspond to N = 0.64 in our 
model) decreasing of the ground state energy due to a 
trapping by the bending with k = 1 (k = 3) comprise 
10% (46%) of the whole excitation energy. This may 
be also important for the dynamics of DNA molecule 
where as was mentioned above the vibration-excitations 
are described by Eq. ([[]) with the dipole-dipole disper- 
sive interaction p{| . Finally, a violation of the Vakhitov- 
Kolokolov stability criterion in the case when the insta- 
bility is due to the softening of the Peierls internal mode 
is found. 
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